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Abstract 

For the first time, differential inclusive-jet cross sections have been measured 
in neutral current deep inelastic ep scattering using the anti-Zcy and SIScone 
algorithms. The measurements were made for boson virtualities Q 2 > 125 GeV 2 
with the ZEUS detector at HERA using an integrated luminosity of 82 pb _1 
and the jets were identified in the Breit frame. The performance and suitability 
of the jet algorithms for their use in hadron-like reactions were investigated by 
comparing the measurements to those performed with the &t algorithm. Next- 
to-leading-order QCD calculations give a good description of the measurements. 
Measurements of the ratios of cross sections using different jet algorithms are 
also presented; the measured ratios are well described by calculations including 
up to 0(a^) terms. Values of a s (Mz) were extracted from the data; the results 
are compatible with and have similar precision to the value extracted from the 
kx analysis. 
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1 Introduction 



Jet production at high transverse energies in neutral-current (NC) deep inelastic ep scat- 
tering (DIS) at HERA provides unique tests of perturbative QCD (pQCD) in a cleaner 
hadron-like enviroment than that encountered in hadron-hadron collisions. In NC DIS, 
the jet search is performed in the Breit frame [T] in which, at leading order (LO) in the 
strong coupling constant, a s , the boson-gluon-fusion (V*g —> qq, with V = 7, Z°) and 
QCD-Compton (V*q — > qg) processes give rise to two hard jets with opposite transverse 
momenta. 

In previous publications, the observables used to test pQCD included inclusive-jet [2]- 
H], dijet jUEHSHU] and multijet [H1II2HII] cross sections and the internal structure of 
jets pT5HT8] defined using the kx cluster algorithm [19] in the longitudinally invariant 
inclusive mode |20j. These studies demonstrated that this jet algorithm results in the 
smallest uncertainties in the reconstruction of jets in ep collisions. The kx algorithm 
is well suited for ep collisions and yields infrared- and collinear-safe cross sections at 
any order of pQCD. However, it might not be best suited to reconstruct jets in hadron- 
hadron collisions, such as those at the LHC. In order to optimise the reconstruction of 
jet observables in such environments, new algorithms were recently developed, like the 
anti-fcr [21], a recombination- type jet algorithm, and the "Seedless Infrared-Safe" cone 
(SIScone) [22] algorithms. The measurements of jet cross sections in NC DIS using these 
algorithms provide additional tests of QCD and also a test of their performance with data 
in a well understood hadron- induced reaction. 

This letter presents measurements of differential inclusive-jet cross sections as a function 
of the jet transverse energy in the Breit frame, E^ B , and the virtuality of the exchanged 
boson, Q 2 , based on the anti-fc^ and SIScone jet algorithms. Differential inclusive-jet cross 
sections as functions of E^ B in different regions of Q 2 are also presented. The analysis is 
based on the same data sample which was used in recent publications [H|5] to make cross- 
section measurements with the kx jet algorithm. The results are compared with next- 
to-leading-order (NLO) QCD calculations using recent parameterisations of the parton 
distribution functions (PDFs) of the proton [23H2B] and with previous measurements 
based on the kx algorithm. Furthermore, measurements of the ratios of inclusive-jet cross 
sections using different jet algorithms are presented and compared to pQCD calculations 
including up to 0(oP s ) terms. In addition, new determinations of a s (Mz) were obtained 
so as to quantify the performance of the anti-^T and SIScone algorithms in comparison 
with that of the kx analysis. 
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2 Jet algorithms 



The study of jet production in hadronic reactions has been well established as a testing 
ground of pQCD. Such tests rely on a trustworthy reconstruction of the topology of the 
final-state partons. This is best done by using jet algorithms with analogous implemen- 
tation in experiment and theory. On the experimental side, the results of the application 
of a jet algorithm should not depend significantly on the presence of soft particles or 
particles which undergo strong decays. On the theoretical side, they should be infrared 
and collinear safe, so that order-by-order pQCD calculations can be performed. A close 
correspondence between the jets and the final-state partons, so that hadronisation cor- 
rections are small, and suppression of beam-remnant jet contributions, when necessary, 
are also required. 

There are two different methods to reconstruct jets from the final-state particles: cluster- 
and cone-type algorithms. Cluster algorithms, such as JADE [26] or kr, are based on 
successive recombinations of particles. They are usually applied in e + e~ experiments, 
where the initial state is governed by QED and the final-state hadrons arise uniquely 
from the short-distance interaction so that all hadrons observed should be associated 
with the hard process and thus clustered. Cone algorithms, such as the iterative cone 
algorithm [27], are based on a maximisation of the energy density within a cone of fixed 
size and are usually applied in hadron collisions, where the initial state consists of coloured 
partons. The initial parton carries only a fraction of the momentum of the parent hadron, 
and the spectator partons lead to the presence of remnant jets and the underlying event, 
which leads to the production of soft hadrons not correlated with the hard interaction. 
Thus, not all final-state particles would be associated to jets. 

In NC DIS, jets are usually defined using the transverse-energy flow in the pseudorapidity 
(?7b)— azimuth (0 B ) plane of the Breit frame. The procedure to reconstruct jets with the 
k,T algorithm from an initial list of objects (e.g. final-state partons, final-state hadrons or 
energy deposits in the calorimeter) is described below in some detail. In the following dis- 
cussion, Ej<b denotes the transverse energy, rf B the pseudorapidity and <p B the azimuthal 
angle of object % in the Breit frame. For each pair of objects, the quantity 

dy = min((^ B ) 2 , (£f B ) 2 ) • [(ri - Vb ) 2 + (0 B - <j> B ) 2 ]/R 2 (1) 

is calculated. For each individual object, the distance to the beam, = (E^ B ) 2 , is 
also calculated. If, of all the values {dij,di}, dki is the smallest, then objects k and I 
are combined into a single new object. If, however, c4 is the smallest, then object k is 
considered a jet and removed from the sample. The procedure is repeated until all objects 
are assigned to jets. 
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The kr algorithm is infrared and collinear safe to all orders in pQCD. As a result of 
the distance between objects defined in Eq. (1), the jets thus obtained have irregular 
shapes [28], in contrast to jets defined using a cone algorithm. From an experimental 
point of view, it may be desirable for jets to have a well-defined shape, e.g. to calibrate the 
jet-energy scale or to estimate the underlying-event contribution. On the other hand, cone 
algorithms produce jets that have an approximate circular shape in the pseudorapidity- 
azimuth plane. However, most implementations of cone algorithms rely on "seeds" to 
start the search for stable cones, which results in unsafe infrared and collinear behaviour 
beyond a given order in pQCD. 

Recently, new jet algorithms have been proposed which produce jets with an approximate 
circular shape but maintain the infrared and collinear safety to all orders in pQCD. Two 
approaches have been followed. The SIScone jet algorithm has been developed to overcome 
the seed problem. A new recombination-type algorithm, the ant i- fey, has been devised 
which, after a modification of the distance as defined in Eq. (1), yields circular jets. 

The SIScone algorithm consists of two steps. First, for a given set of initial objects, all 
stable cones are identified; cones are classified as stable by the coincidence of the cone axis 
with that defined by the total momentum of the objects contained in the given cone of 
radius R in the rj B — B plane of the Breit frame. In this procedure, no seed is used. Stable 
cones are then discarded if their transverse momentum is below a given threshold, pt, m m- 
For each selected stable cone, the scalar sum of the transverse momentum of the objects 
associated to it, p t , is defined. Second, overlapping cones are identified and subsequently 
split or merged according to the following procedure. Two cones are merged if the scalar 
sum of the transverse momentum of the objects shared by the two cones exceeds a certain 
fraction / of the lowest-pt cone; otherwise, two different cones are considered and the 
common objects are assigned to the nearest cone. This jet algorithm is infrared and 
collinear safe to all orders in pQCD. 

The k? and anti-fc^ algorithms are identical except for a modified distance measure, 

dy = min((^ B )- 2 , (^ >B )" 2 ) • [(r? B - r? B ) 2 + (0 B - &]/R 2 , (2) 

and the distance to the beam, which is defined as di = (Ej, B ) 2 - This procedure preserves 
the infrared and collinear safety of the fcr algorithm and, in addition, gives rise to circular 
jets. 

For the measurements presented in this letter, the parameter R was set to unity and the 
jet variables were defined according to the Snowmass convention In the application 
of the SIScone algorithm, the fraction / was set to 0.75 and the transverse momentum 
threshold pt,mm was set to zero. 
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3 Experimental set-up 



A detailed description of the ZEUS detector can be found elsewhere [3U1EI]- A brief 
outline of the components that are most relevant for this analysis is given below. 

Charged particles were tracked in the central tracking detector (CTD) [32], which operated 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consisted of 72 cylindrical drift-chamber layers, organised in nine superlayers covering the 
polar-angld^ region 15° < 9 < 164°. The transverse-momentum resolution for full-length 
tracks can be parameterised as cr(p T )/p T = 0.0058pr © 0.0065 © 0.0014/p T , with p T in 
GeV. The tracking system was used to measure the interaction vertex with a typical 
resolution along (transverse to) the beam direction of 0.4 (0.1) cm and to cross-check the 
energy scale of the calorimeter. 

The high- resolution uranium-scintillator calorimeter (CAL) [33] covered 99.7% of the 
total solid angle and consisted of three parts: the forward (FCAL), the barrel (BCAL) 
and the rear (RCAL) calorimeters. Each part was subdivided transversely into towers and 
longitudinally into one electromagnetic section and either one (in RCAL) or two (in BCAL 
and FCAL) hadronic sections. The smallest subdivision of the calorimeter was called a 
cell. Under test-beam conditions, the CAL single-particle relative energy resolutions were 
a(E)/E = 0.18/ VE for electrons and a(E)/E = 0.35 /VE for hadrons, with E in GeV. 

The luminosity was measured from the rate of the bremsstrahlung process ep — > e'jp. The 
resulting small-angle energetic photons were measured by the luminosity monitor [31], a 
lead-scintillator calorimeter placed in the HERA tunnel at Z = — 107 m. 



4 Data selection 

The data were collected during the running period 1998-2000, when HERA operated with 
protons of energy E p = 920 GeV and electron^] of energy E e = 27.5 GeV, and correspond 
to an integrated luminosity of 81. 7± 1.8 pb" 1 , of which 16.7 pb^ 1 (65.0 pb _1 ) was for e~p 
(e + p) collisions. 

Neutral current DIS events were selected using the same criteria as reported in recent 
publications [HE]. The main steps are briefly listed below. 

1 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing towards the 
centre of HERA. The coordinate origin is at the nominal interaction point. 

2 Here and in the following, the term "electron" denotes generically both the electron (e~) and the 
positron (e + ). 
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The scattered-electron candidate was identified from the pattern of energy deposits in the 
CAL [53]. The energy (E' e ) and polar angle of the electron candidate were determined 
from the CAL measurements. The Q 2 variable was reconstructed using the double-angle 
method [36]. The angle 7^, which corresponds to the angle of the scattered quark in the 
quark-parton model, was reconstructed using the hadronic final state [36] . 

The main requirements imposed on the data sample were: an electron candidate with 
E' e > 10 GeV; a vertex position along the beam axis in the range \Z\ < 34 cm; 38 < 
(E — Pz) < 65 GeV, where E is the total energy as measured by the CAL, E = 'Yl li E i , 
and Pz is the Z component of the vector P = Eiff, in both cases the sum runs over 
all CAL cells, Ei is the energy of the CAL cell % and f, is a unit vector along the line 
joining the reconstructed vertex and the geometric centre of the cell i; Q 2 > 125 GeV 2 ; 
and I cos 7/^ < 0.65. In this selected sample, contamination from non-ep interactions and 
other physics processes was negligible. 

The anti-A/r and SIScone algorithms^ were both used to reconstruct jets in the hadronic 
final state both in data and in Monte Carlo (MC) simulated events (see Section EJ. In the 
data, the algorithms were applied to the energy deposits measured in the CAL cells after 
excluding those associated with the scattered-electron candidate. After reconstructing 
the jet variables in the Breit frame, the massless four-momenta were boosted into the 
laboratory frame, where the transverse energy (-E^lab) an d the pseudorapidity (^lab) °f 
each jet were calculated. Energy corrections j2j[T5j|38] were then applied to the jets in 
the laboratory frame, separately for the anti-fc^ and SIScone jet samples. In addition, 
events were removed from the samples if any of the jets was in the backward region of the 
detector (r}^ AB < —2) and jets were not included in the final samples if E^ LAB < 2.5 GeV. 

Only events with at least one jet in the pseudorapidity range —2 < rj B < 1.5 were kept 
for further analysis. The final data samples with at least one jet satisfying E^ B > 8 GeV 
contained 18847 events (12575 one-jet, 6126 two-jet, 145 three-jet and one four-jet) in 
the anti-A/r sample and 19486 events (13081 one-jet, 6252 two-jet, 152 three-jet and one 
four-jet) in the SIScone sample. 



5 Monte Carlo simulation 

Samples of events were generated to determine the response of the detector to jets of 
hadrons and to compute the correction factors necessary to obtain the hadron-level jet 
cross sections. The hadron level is defined by those hadrons with lifetime r > 10 ps. 

3 The implementations of the anti-fcr and SIScone algorithms in the Fastjet package (37], version 2.4.1, 
were used. 
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The generated events were passed through the Geant 3.13-based (39] ZEUS detector- 
and trigger-simulation programs [31]. They were reconstructed and analysed by the same 
program chain as the data. 

Neutral current DIS events including electroweak radiative effects were simulated using 
the Heracles 4.6.1 [3D] program with the Djangoh 1.1 [31] interface to the QCD 
programs. The QCD cascade was simulated using the colour-dipole model (CDM) [32] 
including the LO QCD diagrams as implemented in Ariadne 4.08 [33] and, alternatively, 
with the MEPS model of Lepto 6.5 [44] . The CTEQ5D [3S] parameterisations of the 
proton PDFs were used for these simulations. Fragmentation into hadrons was performed 
using the Lund string model [3D] as implemented in Jetset [4TII45] . 

The jet search was performed on the simulated events using the energy measured in the 
CAL cells in the same way as for the data. The same jet algorithms were also applied 
to the final-state particles (hadron level) and to the partons available after the parton 
shower (parton level). Additional MC samples were used to correct the measured cross 
sections for QED radiative effects and the running of a cm . 

6 Acceptance corrections and experimental uncer- 
tainties 

To measure the cross sections, the E^ B and Q 2 distributions in the data were corrected 
for detector effects using bin-by-bin correction factors determined with the MC samples. 
These correction factors took into account the efficiency of the selection criteria and 
the purity and efficiency of the jet reconstruction. For this approach to be valid, the 
uncorrected distributions of the data must be well described by the MC simulations at the 
detector level. This condition was satisfied by both the Ariadne and Lepto-MEPS MC 
samples. The average between the acceptance-correction values obtained with Ariadne 
and Lepto-MEPS was used to correct the data to the hadron level. The deviations in 
the results obtained by using either Ariadne or Lepto-MEPS to correct the data from 
their average were taken to represent systematic uncertainties of the effect of the QCD- 
cascade model in the corrections (see below). The acceptance-correction factors differed 
from unity by typically less than 10%. 

The following sources of systematic uncertainty were considered for the measured cross 
sections: 

• the uncertainty in the absolute energy scale of the jets was estimated to be ±1% 
for £^ L AB > 10 GeV and ±3% for lower E^ LAB values [381 [Ml [50] - The resulting 
uncertainty on the cross sections was about ±5%; 



6 



• the differences in the results obtained by using either Ariadne or Lepto-MEPS to 
correct the data for detector effects were taken to represent systematic uncertainties. 
The resulting uncertainty on the cross sections was typically below ±3%; 

• the uncertainty due to the selection cuts was estimated by varying the values of the 
cuts within the resolution of each variable; the effect on the cross sections was typically 
below ±3%; 

• the uncertainty on the reconstruction of the boost to the Breit frame was estimated 
by using the direction of the track associated with the scattered electron instead of 
that derived from its impact position in the CAL. The effect on the cross sections was 
typically below ±1%; 

• the uncertainty in the absolute energy scale of the electron candidate was estimated 
to be ±1% [51]. The resulting uncertainty on the cross sections was below ±1%; 

• the uncertainty in the cross sections due to that in the simulation of the trigger was 
negligible. 

The systematic uncertainties were similar for both jet algorithms. Those not associated 
with the absolute energy scale of the jets were added in quadrature to the statistical 
uncertainties and are shown in the figures as error bars. The uncertainty due to the 
absolute energy scale of the jets is shown separately as a shaded band due to the large 
bin-to-bin correlation. In addition, there was an overall normalisation uncertainty of 2.2% 
from the luminosity determination, which is not included in the figures. 



7 QCD calculations and theoretical uncertainties 

Next-to-leading-order (0(a 2 s )) QCD calculations were compared to the measured inclusive- 
jet cross sections. The calculations were obtained using the program Disent [52], in the 
MS renormalisation and factorisation schemes using a generalised version [52] of the sub- 
traction method [53]. The number of flavours was set to five and the renormalisation (//r) 
and factorisation (jip) scales were chosen to be [Mr = E^ B and pip — Q, respectively. The 
strong coupling constant was calculated at two loops with = 226 MeV, corresponding 
to a s (M z ) = 0.118. The calculations were performed using the ZEUS-S [23] parameter- 
isations of the proton PDFs^l In Disent, the value of a cm was fixed to 1/137. The 
anti-fcr and SIScone algorithms were also applied to the partons in the events generated 
by Disent in order to compute the jet cross-section predictions. 

4 The LHAPDF package [S3], version 5.7.1, was used. 
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Measurements are also presented of the ratio of inclusive-jet cross sections based on differ- 
ent jet algorithms. Although the calculations for inclusive-jet cross sections can be made 
at present only up to 0(a 2 ) (see Fig. [TJ, the difference^ of cross sections using different 
jet algorithms can be predicted up to 0(a z s ) using the program NlojetH — \- [55]. The 



same parameter settings as defined above were used for the Nlojet-I — \- predictions. The 
predicted ratio (da Sknt i-k T /dX)/(dak T /dX), where X = Q 2 or E^ B , was calculated as 



dCanti-k T /dX _ _|_ d<J anti _ kT / dX — d<Jk T /dX ^ ^ _|_ C ' a l /g\ 
da kT /dX da kr /dX A ■ a s + B ■ a 2 ' 

since differences between the anti-fey and fey algorithms appear first for final states with 
four partons and were evaluated using the tree- level eq — > eqggg, eq —> eqgqq, eg — > eqqqq 
and eg — > eggqq subprocesses (see Fig.[T)). The predicted ratio (dasiscone/ dX) / (da kT / dX) 
was calculated as 



dasisconc/dX dasiscone / dX - da kT /dX D ■ a 2 s + E ■ a 3 s 

= 1 T — 1 ~T~ (4) 

dak T /dX dak T / dX A ■ a s + B ■ a 2 ' 

since differences between the SIScone and kr algorithms appear first for final states with 
three partons and were evaluated using the tree-level three-parton and four-parton sub- 
processes and the one-loop three-parton configurations (see Fig. [T]). Equation (4) also 
applies to the predicted ratio (da an ti-k T / dX) / (dasiscone/ dX) . 

Since the measurements refer to jets of hadrons, whereas the NLO QCD calculations refer 
to jets of partons, the predictions were corrected to the hadron level using the MC models. 
The multiplicative correction factor (Chad) was defined as the ratio of the cross section for 
jets of hadrons over that for jets of partons, estimated by using the MC programs described 
in Section [51 The mean of the ratios obtained with Ariadne and Lepto-MEPS was 
taken as the value of Chad- The value of Chad for the inclusive-jet cross sections differs 
from unity by less than 5% [5], 6% and 11% for the kr, anti-fey and SIScone jet algorithms, 
respectively, in the region Q 2 > 500 GeV 2 . 

Neither Disent nor NlojetH — \- include the contribution from Z° exchange; MC simu- 
lated events with and without Z° exchange were used to include this effect in the pQCD 
predictions. In the following, pQCD calculations refer to the fully corrected predictions. 

Several sources of uncertainty in the theoretical predictions for the inclusive-jet cross 
sections were considered: 

• the uncertainty on the NLO QCD calculations due to terms beyond NLO, estimated 
by varying /!# between E^ B /2 and 2£ T it B , was below ±7 (±10)% at low Q 2 and low 

5 The differences of cross sections were evaluated on an event- by- event basis. 
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E^ B and decreased to less than ±5 (±7)% for Q 2 > 250 GeV 2 for the k T [5] and 
anti-kx (SIScone) algorithms; 

• the uncertainty on the NLO QCD calculations due to those on the proton PDFs was 
estimated by repeating the calculations using 22 additional sets from the ZEUS-S anal- 
ysis, which takes into account the statistical and correlated systematic experimental 
uncertainties of each data set used in the determination of the proton PDFs. The 
resulting uncertainty in the cross sections for all three jet algorithms was below ±3%, 
except in the high-E^g region where it reached ±4%; 

• the uncertainty on the NLO QCD calculations due to that on a s {Mz) was estimated by 
repeating the calculations using two additional sets of proton PDFs, for which different 
values of a s [Mz) were assumed in the fits. The difference between the calculations 
using these various sets was scaled by a factor to reflect the uncertainty on the current 
world average of a s [SS]- The resulting uncertainty in the cross sections was below 
±2% for all three jet algorithms; 

• the uncertainty from the modelling of the QCD cascade was assumed to be half the 
difference between the hadronisation corrections obtained using the Ariadne and 
Lepto-MEPS models. The resulting uncertainty on the cross sections was less than 
±1.4% [5], ±1.7% and ±2.3% for the kr, anti-£>r and SIScone algorithms, respectively; 

• the uncertainty of the calculations due to the choice of /ip was estimated by repeating 
the calculations with fip = Q/2 and 2Q. The effect was negligible. 

The total theoretical uncertainty was obtained by adding in quadrature the individual 
uncertainties listed above. As a function of Q 2 , the total theoretical uncertainty varies in 
the range 3 — 7% (3 — 10%) for the anti-fc T (SIScone) algorithm; as a function of -E^g, 
the range of variation is 5 — 6% (6 — 8%) for the anti-£>r (SIScone) algorithm. 

It is concluded that the NLO QCD predictions for the inclusive-jet cross sections are 
of similar precision for the kx and anti-fcy algorithms and somewhat less precise for the 
SIScone. A more detailed comparison of the theoretical uncertainties is presented in 
Section |HJ 

8 Results 

8.1 Inclusive-jet differential cross sections with different jet al- 
gorithms 

The inclusive-jet differential cross sections were measured in the kinematic region Q 2 > 
125 GeV 2 and |cos7/j| < 0.65. The jets were reconstructed using either the anti-A/p or 
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the SIScone jet algorithms and the cross sections refer to jets with E 3 ^ B > 8 GeV and 
—2 < rf B < 1.5. These cross sections were corrected for detector and QED radiative 
effects and the running of a em . 

The measurements of the inclusive-jet differential cross sections as functions of E^ B and 
Q 2 are presented in Fig. [2] and Tables [T] and [2] for the anti-fey and SIScone jet algo- 
rithms. For comparison, the measurements of inclusive-jet cross sections with the fey 
algorithm [5] are also shown. In Fig. El each data point is plotted at the abscissa at which 
the NLO QCD differential cross section was equal to its bin-averaged value. The measured 
da/dEjf-g (da/dQ 2 ) exhibits a steep fall-off over three (five) orders of magnitude for the 
jet algorithms considered in the E^ B (Q 2 ) measured range. The measurements using the 
three jet algorithms have a very similar shape and normalisation. 

The NLO QCD predictions with fiR = E^ B are compared to the measurements in Fig. [2j 
The hadronisation-correction factors and their uncertainties are also shown. The ratios 
of the measured differential cross sections to the NLO QCD calculations are shown in 
Figs. Eh and [3b separately for each jet algorithm. The measured differential cross sections 
are well reproduced by the calculations, with similar precision in all cases. To study the 
scale dependence, NLO QCD calculations using fiR = Q were also compared to the data 
(not shown); they also provide a good description of the data. 

To compare in more detail the results of the different algorithms, the ratios of the measured 
cross-sections anti-fey/fey, SIScone/ fey and anti- fey /SIScone were investigated (see Figs. 
and[3H). In the ratios, the statistical correlations among the event samples as well as 
those among the jets in the same event were taken into account in the estimation of 
the statistical uncertainties. These ratios show that the cross sections with the three jet 
algorithms are similar: the ratios as functions of Q 2 differ from unity by less than 3.2%; as 
functions of E^g, they differ by less than 3.6% except at high S^g, where the difference 
is 10%. The pQCD predictions including up to 0{o? s ) terms for the ratios of cross sections 
are also shown in Figs. &p and Eli. In the estimation of the total theoretical uncertainty 
of the predicted ratio, all the theoretical contributions were assumed to be correlated 
except that due to terms beyond 0{pt z s ) and that from the modelling of the QCD cascade. 
The measured ratios are well described by the calculations including up to O(af) terms 
within the small experimental and theoretical uncertainties, which are dominated by the 
uncertainty on the modelling of the QCD cascade. 

Figure H] shows the contributions to the theoretical uncertainty of the inclusive-jet cross 
sections from terms beyond NLO, from the uncertainty in the PDFs, from that on the 
value of a s (Mz) and from that on the modelling of the QCD cascade, separately for the 
anti-fey and SIScone algorithms. The uncertainty coming from the terms beyond NLO is 
dominant at low -Ey B and low Q 2 in all cases and somewhat higher for the predictions 
based on SIScone than those for the anti-fey and fey [5] algorithms. At high E^ B , the PDF 
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uncertainty is dominant in the case of the anti-fcy algorithm and of the same order as that 
arising from terms beyond NLO for the SIScone algorithm. At high Q 2 , the uncertainty 
due to terms beyond NLO is of the same order for both algorithms and remains dominant. 
The uncertainty on the modelling of the QCD cascade is somewhat higher for the SIScone 
algorithm than for anti-A>r, but significantly smaller than that from terms beyond NLO. 
The total theoretical uncertainty for the predictions using //r = Q is slightly larger than 
that using fi R = E 3 ^ B . 

The inclusive-jet cross sections for different regions of Q 2 as a function of -E^g are pre- 
sented in Fig. |5] and Tables |3] and H] for both algorithms. The measured cross sections 
exhibit a steep fall-off within the -Ej^ B range considered. As Q 2 increases, the E 3 ^ B de- 
pendence of the cross section becomes less steep. The measurements have a similar shape 
and normalisation for both jet algorithms and are similar to the results obtained with 
the algorithm [1]. The NLO QCD predictions are compared to the measurements in 
Fig. |5j The data are well described by the predictions. Figure |6] shows the ratios of the 
measured differential cross sections to the NLO QCD calculations. The uncertainty of 
the NLO QCD calculations is also shown: the uncertainties from the anti-fcr predictions 
are of a similar size as those encountered for the kx algorithm |3], and those for SIScone 
are somewhat larger at low E 3 ^ and low Q 2 . 

In summary, it is concluded that the data for the fcy, anti-fcy and SIScone jet algorithms 
are well described by the NLO QCD calculations with similar experimental and theoret- 
ical precision. Furthermore, the measured ratios are well described by the predictions 
including up to 0(oP s ) terms, demonstrating the ability of the pQCD calculations with 
up to four partons in the final state to account adequately for the differences between the 
jet algorithms. 

8.2 Determination of ot s (Mz) 

The measured differential cross-sections da/dQ 2 for Q 2 > 500 GeV 2 were used to de- 
termine values of a s {M z ) using the method presented previously [2]. The NLO QCD 
calculations were performed using the program Disent with five sets of ZEUS-S proton 
PDFs which were determined from global fits assuming different values of a s (Mz), namely 
a s {M z ) = 0.115, 0.117, 0.119, 0.121 and 0.123. The value of a s {M z ) used in each calcu- 
lation was that associated with the corresponding set of PDFs. The a s (Mz) dependence 
of the predicted cross sections in each bin i of Q 2 was parameterised according to 

[da/dQ 2 (a s (M z ))] i = C{a s {M z ) + C\a 2 s {M z ), 
where C\ and C % 2 were determined from a \ 2 fit to the NLO QCD calculations. The value 
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of a s [Mz) was determined by a x 2 fit to the measured da/dQ 2 values. 

The uncertainties on the extracted values of a s (Mz) due to the experimental systematic 
uncertainties were evaluated by repeating the analysis for each systematic check presented 
in Section [6j The overall normalisation uncertainty from the luminosity determination 
was also considered. The largest contribution to the experimental uncertainty comes 
from the jet energy scale and amounts to ±1.9% on a s (Mz) for both algorithms (±2% for 
fcr U). The theoretical uncertainties were evaluated as described in Section[7J The largest 
contribution was the theoretical uncertainty arising from terms beyond NLO, which was 
estimated by using the method proposed by Jones et al. [57], and amounted to ±1.5% 
for both algorithm^] (±1.5% for kx [5]). The uncertainty due to the proton PDFs was 
±0.7 (±0.8)% for the anti-fcr (SIScone) algorithm (±0.7% for kr [5]). The uncertainty 
arising from the hadronisation effects amounted to ±0.9 (±1.2)% for the anti-Zc^ (SIScone) 
algorithm (±0.8% for kx [5]). Thus, the performance of the three jet algorithms is very 
similar. 

As a cross-check, a s [Mz) was determined by using NLO QCD calculations based on the 
CTEQ6.1 [24] (MSTW2008 [25]) sets of proton PDFs: the values obtained are consistent 
within ±1.2 (±1.0)% with those based on ZEUS-S. The uncertainty arising from the 
proton PDFs was estimated to be ±1.5 (±0.7, ±0.4)% using the results of the CTEQ6.1 
(MSTW2008nlo90cl, MSTW2008nlo68cl) analysis. 

The values of a s (Mz) obtained from the measured da/dQ 2 for Q 2 > 500 GeV 2 are 



«s(M z )| anti _ fcT = 0.1188 ±0.0014 (stat.) l^ol (exp.) Ijffi (th.) and 
« s (M z )|siscone = 0.1186 ± 0.0013 (stat.) (exp.) Ijffi (th.) . 

These values of a s (Mz) are consistent with that obtained using the kx algorithm [5], 

a s (M z )\ kT = 0.1207 ±0.0014 (stat.) t°;°°j* (exp.) (th.), 

with the current world average of 0.1189 ± 0.0010 [56], with the results obtained in pp 
collisions [58] as well as with the HERA average of 0.1186±0.0051 [59]. It should be noted 
that the differences between the central values of a s (Mz) obtained using the three jet algo- 
rithms are comparable to the uncertainties due to higher-order terms in the calculations. 
It is observed that the precision in a s [Mz) obtained with the k?, anti-fcy and SIScone jet 
algorithms is very similar, and comparable to those obtained in e + e~ interactions [56J. 

6 The theoretical uncertainty arising from terms beyond NLO estimated by refitting the data using 
calculations based on = 2JByg or E^ t B /2 amounts to Ig'J (ig'y)% for the anti-fcr (SIScone) 
algorithm. 
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9 Summary and conclusions 



For the first time, differential cross sections for inclusive-jet production in neutral current 
deep inelastic ep scattering were measured using the anti-fey and SIScone jet algorithms 
with R — 1. The cross sections correspond to a centre-of-mass energy of 318 GeV and 
refer to jets of hadrons with > 8 GeV and —2 < r) B < 1.5 identified in the Breit 
frame with the anti-&x or SIScone jet algorithms. The measurements are given in the 
kinematic region of Q 2 > 125 GeV 2 and | cos 77^ < 0.65. 

A detailed comparison between the measurements as functions of E^ B and Q 2 for both 
algorithms and those from a previous analysis based on the kx jet reconstruction was 
performed. The measured cross sections for the three jet algorithms have similar shapes 
and normalisations. 

The NLO QCD calculations of inclusive-jet cross sections and their uncertainties for the 
different algorithms were also compared: the data are well described by the predictions; 
the calculations for the anti-fey algorithm have a similar precision as that of the kx whereas 
those for the SIScone are somewhat less precise due to the contribution from terms beyond 
NLO. 

Measurements of the ratios of cross sections using different jet algorithms were also pre- 
sented and compared to calculations including up to 0(a%) terms. The measured ratios 
are well reproduced by the predictions, demonstrating the ability of the pQCD calculations 
including up to four partons in the final state to account adequately for the differences 
between the jet algorithms. 

Values of a s (Mz) were extracted from the measured inclusive-jet differential cross sections 
using each algorithm. QCD fits of the cross-sections da/dQ 2 for Q 2 > 500 GeV 2 yield the 
following values of a s (Mz): 



a s (M z )| anti _ fcT = 0.1188 ±0.0014 (stat.) (exp.) Iffi (th.) and 

a s (M z )| SISconc = 0.U86 ±0.0013 (stat.) (exp.) i^ol (th-). 

These values are consistent with each other and with that obtained from the kx analysis 
with a similar precision. 
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Ej, B bin 


da 1 dE^ B 


"stat 


^syst 


<5es 


Cqed 


Chad • Czo 


(GeV) 


(pb/GeV) 













anti-&T 



8- 


- 10 


62.84 


±0.69 


+1.40 
-1.49 


+2.67 
-2.46 


0.95 


0.93 


10 


- 14 


28.56 


±0.34 


+0.39 
-0.48 


+1.31 
-1.25 


0.97 


0.93 


14 


- 18 


10.69 


±0.20 


+0.20 
-0.23 


+0.57 
-0.49 


0.96 


0.92 


18 


-25 


3.186 


±0.082 


+0.041 
-0.040 


+0.143 
-0.154 


0.95 


0.92 


25 


-35 


0.725 


±0.032 


+0.014 
-0.014 


+0.038 
-0.034 


0.94 


0.94 


35- 


- 100 


0.0300 


±0.0026 


+0.0005 
-0.0006 


+0.0014 
-0.0015 


1.06 


0.92 



SIScone 



8- 


- 10 


64.08 


±0.69 


+1.74 
-1.82 


+2.73 
-2.54 


0.95 


0. 


B7 


10 


- 14 


28.92 


±0.34 


+0.42 
-0.49 


+1.42 
-1.29 


0.96 


0. 


35 


14 


- 18 


10.72 


±0.20 


+0.16 
-0.22 


+0.55 
-0.51 


0.96 


0. 


B5 


18 


-25 


3.211 


±0.081 


+0.045 
-0.043 


+0.161 
-0.149 


0.95 


0. 


B5 


25 


-35 


0.744 


±0.033 


+0.017 
-0.017 


+0.036 
-0.039 


0.93 


0. 


B8 


35- 


- 100 


0.0332 


±0.0028 


+0.0008 
-0.0008 


+0.0020 
-0.0016 


1.06 


0. 


B6 



Table 1: The measured differential cross-sections da/dE^ B for inclusive-jet 
production using different jet algorithms. The statistical, uncorrelated systematic 
and jet-energy-scale (ES) uncertainties are shown separately. The multiplicative 
corrections applied to the data to correct for QED radiative effects, Cqed ; and the 
corrections for hadronisation and Z° effects to be applied to the parton-level NLO 
QCD calculations, C*h a d ■ Czo, are shown in the last two columns. 
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Q 2 bin 


da/dQ 2 


"stat 


<^syst 


<$ES 


Cqed 


Chad ' C z o 


(GeV 2 ) 


(pb/GeV 2 ) 













anti-fc^ 



125 


- 250 


1.078 


±0.012 


+0.017 
-0.018 


+0.067 
-0.060 


0.97 


0.90 


250 


-500 


0.3589 


±0.0053 


+0.0055 
-0.0084 


+0.0147 
-0.0142 


0.95 


0.94 


500 


- 1000 


0.1004 


±0.0020 


+0.0019 
-0.0019 


+0.0030 
-0.0028 


0.95 


0.94 


1000 


- 2000 


0.02406 


±0.00072 


+0.00044 
-0.00033 


+0.00045 
-0.00058 


0.94 


0.96 


2000 


- 5000 


0.00385 


±0.00017 


+0.00015 
-0.00012 


+0.00007 
-0.00007 


0.94 


0.97 


5000 - 


- 100000 


36.6 ■ 10^ 6 


±3.3 ■ 10~ 6 


+2.2 1n -6 
-2.6 ' iu 


+ 1.0 . ln -6 
-0.9 1U 


0.98 


0.95 



SIScone 



125 


-250 


1.099 


±0.012 


+0.020 
-0.021 


+0.069 
-0.063 


0.97 


0.82 


250 


-500 


0.3637 


±0.0053 


+0.0062 
-0.0089 


+0.0156 
-0.0143 


0.95 


0.87 


500 


- 1000 


0.1011 


±0.0020 


+0.0025 
-0.0024 


+0.0030 
-0.0028 


0.95 


0.89 


1000 


- 2000 


0.02419 


±0.00071 


+0.00048 
-0.00045 


+0.00060 
-0.00057 


0.94 


0.92 


2000 


- 5000 


0.00393 


±0.00017 


+0.00014 
-0.00012 


+0.00007 
-0.00007 


0.94 


0.93 


5000 - 


- 100000 


36.8 ■ 10^ 6 


±3.2 • 10~ 6 


l£g • io- 6 


+ 1.1 . in-6 
-0.6 1U 


0.97 


0.92 



Table 2: The measured differential cross-sections da/dQ 2 for inclusive-jet pro- 
duction using different jet algorithms. Other details as in the caption to TableUl 
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4% 


bin 


7 / 7 7 — ll Ct 

da j ah J T B 


fstat 


Csyst 




<^QED 


L'had ■ <^Z° 


(GeV) 


(pb/GeV) 


















125 < Q 2 


< 250 GeV 2 






8 - 


10 


31.60 


±0.49 


_i_n on 
-0.93 


4-1 R7 
Ti.oi 

-1.56 


0.95 


0.91 


10 - 


14 


12.67 


±0.22 


+0.13 
-0.15 


+0.80 
-0.74 


0.99 


0.90 


14 - 


18 


3.56 


±0.11 


+0.03 
-0.05 


+0.26 
-0.24 


0.96 


0.88 


18 - 


25 


0.790 


±0.037 


+0.013 
-0.012 


+0.051 
-0.047 


0.95 


0.86 


25 — 


100 


0.0165 


±0.0015 


+0.0013 
-0.0013 


+0.0008 
-0.0011 


1.02 


0.85 








250 < Q 2 


< 500 GeV 2 






8 - 


10 


17.90 


±0.38 


+0.41 
-0.50 


+0.52 
-0.58 


0.94 


0.94 


10 - 


14 


8.37 


±0.19 


+0.18 
-0.26 


+0.36 
-0.32 


0.96 


0.94 


14 - 


18 


3.17 


±0.11 


+0.04 
-0.08 


+0.18 
-0.14 


0.97 


0.94 


18 - 


25 


0.882 


±0.042 


+0.016 
-0.015 


+0.047 
-0.055 


0.92 


0.89 


25 - 


100 


0.0225 


±0.0020 


+0.0006 
-0.0006 


+0.0011 
-0.0011 


0.94 


0.95 








500 < Q 


< 1000 GeV 2 






8 - 


10 


8.36 


±0.26 


+0.17 
-0.18 


+0.19 
-0.14 


0.96 


0.95 


10 - 


14 


4.52 


±0.14 


+0.14 
-0.13 


+0.12 
-0.13 


0.94 


0.92 


14 - 


18 


2.116 


±0.092 


+0.099 
-0.097 


+0.069 
-0.058 


0.93 


0.95 


18 - 


25 


0.682 


±0.039 


+0.017 
-0.016 


+0.026 
-0.029 


0.97 


0.95 


25 - 


100 


0.0293 


i r\ r\r\c\ p" 

±0.0025 


+0.0014 
-0.0014 


+0.0019 
-0.0016 


0.94 


0.98 








■i r\r\r\ , / — '\ 

1000 < Q z 


< 2000 GeV 2 






8 - 


10 


3.34 


±0.16 


+0.10 
-0.08 


+0.00 
-0.10 


0.93 


0.97 


10 - 


14 


1.967 


±0.093 


+0.070 
-0.065 


+0.025 
-0.039 


0.93 


0.96 


14 - 


18 


1.099 


±0.068 


+0.049 
-0.048 


+0.036 
-0.015 


0.95 


0.94 


18 - 


25 


0.477 


±0.033 


+0.024 
-0.024 


+0.009 
-0.011 


0.94 


1.00 


25 - 


100 


0.0237 


±0.0023 


+0.0015 
-0.0015 


+0.0015 
-0.0012 


0.99 


0.99 








c\r\r\r\ /~\ 

2000 < Q z 


< 5000 GeV 2 






8 - 


10 


1.374 


±0.099 


+0.094 
-0.084 


+0.034 
-0.004 


0.91 


0.89 


10 - 


14 


0.797 


±0.059 


+0.030 
-0.022 


+0.009 
-0.009 


0.93 


0.96 


14- 


18 


0.598 


±0.051 


+0.032 
-0.026 


+0.025 
-0.012 


1 Ml 

1.01 


0.99 


18 - 


25 


0.261 


±0.026 


+0.006 
-0.012 


+0.005 
-0.007 


0.97 


1.05 


25- 


100 


0.0189 


±0.0022 


+0.0010 
-0.0011 


+0.0007 
-0.0006 


0.90 


1.02 


5000 < Q 2 < 100000 GeV 2 


8- 


10 


0.301 


±0.051 


+0.033 
-0.037 


+0.005 
-0.021 


0.94 


0.89 


10 - 


14 


0.229 


±0.034 


+0.007 
-0.013 


+0.013 
-0.005 


0.96 


0.93 


14- 


18 


0.148 


±0.028 


+0.020 
-0.017 


+0.000 
-0.009 


1.01 


0.92 


18 - 


25 


0.081 


±0.016 


+0.013 
-0.014 


+0.002 
-0.001 


0.98 


0.99 


25- 


100 


0.0104 


±0.0020 


+0.0007 
-0.0007 


+0.0002 
-0.0001 


1.00 


1.04 



Table 3: The measured differential cross-section da/dE 1 ^ for inclusive-jet pro- 
duction in different regions of Q 2 using the anti-kx jet algorithm. Other details as 
in the caption to Tabled 
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4% 


bin 


7 / 7 7 — llCt 

aa/ah J TB 


c - 

"stat 


^syst 


C>ES 


^QED 


L'had ■ <^Z° 


(GeV) 


(pb/GeV) 


















125 < Q 


< 250 GeV 2 






8 - 


10 


32.36 


±0.49 


_i_n qq 
-1.02 


-1-1 ^fi 

-1.67 


0.96 


0.84 


10 - 


14 


12.79 


±0.22 


+0.17 
-0.18 


+0.84 
-0.78 


0.98 


0.81 


14 - 


18 


3.66 


±0.11 


+0.04 
-0.05 


+0.26 
-0.24 


0.97 


0.80 


18 - 


25 


0.792 


±0.037 


+0.009 
-0.008 


+0.056 
-0.044 


0.95 


0.79 


25 — 


100 


0.0180 


±0.0015 


+0.0013 
-0.0013 


+0.0009 
-0.0012 


1.00 


0.80 








250 < Q 2 


< 500 GeV 2 






8 - 


10 


18. 0z 


__i_n ot 

±0.37 


+0.53 
-0.60 


+0.59 
-0.55 


0.94 


0.89 


10 - 


14 


8.55 


±0.19 


+0.19 
-0.24 


+0.37 
-0.33 


0.96 


0.87 


14 - 


18 


3.16 


±0.11 


+0.04 
-0.09 


+0.17 
-0.15 


0.97 


0.86 


18 - 


25 


0.893 


±0.042 


+0.021 
-0.021 


+0.055 
-0.050 


0.92 


0.82 


25 - 


100 


r\ r\c\ A T 

0.0247 


i r\ nnoi 

±0.0021 


+0.0004 
-0.0004 


+0.0012 
-0.0014 


0.94 


0.87 








500 < Q 


< 1000 GeV 2 






8 - 


10 


8.65 


±0.26 


+0.25 
-0.24 


+0.14 
-0.16 


0.96 


0.89 


10 - 


14 


4.49 


±0.14 


+0.12 
-0.11 


+0.13 
-0.12 


0.94 


0.88 


14 - 


18 


2.112 


±0.091 


+0.074 
-0.073 


+0.073 
-0.067 


0.94 


0.90 


18 - 


25 


0.677 


±0.038 


+0.024 
-0.024 


+0.031 
-0.026 


0.97 


0.89 


25 - 


100 


0.0289 


±0.0025 


+0.0010 
-0.0010 


+0.0016 
-0.0018 


0.93 


0.89 








1 r\r\r\ * 

1000 < Q 


< 2000 GeV 2 






8 - 


10 


3.30 


±0.16 


+0.12 
-0.11 


+0.07 
-0.06 


0.92 


0.94 


10 - 


14 


2.024 


±0.096 


+0.080 
-0.085 


+0.038 
-0.036 


0.95 


0.92 


14 - 


18 


1.070 


±0.068 


+0.049 
-0.048 


+0.024 
-0.025 


0.96 


0.91 


18 - 


25 


0.478 


±0.033 


+0.026 
-0.025 


+0.010 
-0.014 


0.93 


0.93 


25 — 


100 


0.0246 


±0.0024 


+0.0018 
-0.0018 


+0.0017 
-0.0012 


0.98 


0.93 








2000 < Q 


< 5000 GeV 2 






8 - 


10 


1.47 


±0.11 


+0.07 
-0.05 


+0.00 
-0.04 


0.94 


0.86 


10 - 


14 


0.822 


±0.059 


+0.037 
-0.031 


+0.020 
-0.010 


0.94 


0.94 


14- 


18 


0.563 


±0.049 


+0.032 
-0.034 


+0.012 
-0.006 


0.97 


0.92 


18 - 


25 


0.276 


±0.027 


+0.010 
-0.015 


+0.006 
-0.012 


0.98 


0.99 


25 - 


100 


0.0186 


±0.0021 


+0.0008 
-0.0008 


+0.0006 
-0.0005 


0.88 


0.98 


5000 < Q 2 < 100000 GeV 2 


8- 


10 


0.264 


±0.044 


+0.047 
-0.052 


+0.014 
-0.001 


0.94 


0.90 


10 - 


14 


0.248 


±0.036 


+0.019 
-0.024 


+0.011 
-0.006 


0.96 


0.90 


14 - 


18 


0.146 


±0.027 


+0.011 
-0.012 


+0.003 
-0.002 


0.99 


0.91 


18 - 


25 


0.082 


±0.016 


+0.015 
-0.016 


+0.002 
-0.002 


1.00 


0.92 


25 - 


100 


0.0111 


±0.0021 


+0.0009 
-0.0008 


+0.0004 
-0.0004 


1.00 


0.99 



Table 4: The measured differential cross-section da/dE 1 ^ for inclusive-jet pro- 
duction in different regions of Q 2 using the SIScone jet algorithm. Other details as 
in the caption to Tabled 
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FINAL STATE PARTONS 



FINAL STATE PARTONS 



FINAL STATE PARTONS 



tree 



Figure 1: Examples of Feynman diagrams contributing to inclusive- jet produc- 
tion in the Breit frame up to 0(a^). The calculations of inclusive- jet cross sections 
up to 0{a 2 s ) include the (a) leading-order diagrams and (b) virtual and (d) real 
corrections. The lowest-order diagrams that contribute to the cross-section differ- 
ence between the anti-kr and kr algorithms are of type (f). The calculations of the 
cross-section difference between the SIScone and kr algorithms include diagrams of 
type (d), (e) and (f). 
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E j x' B > 8 GeV 
-2 < rt>f < 1.5 
Icos y h l < 0.65 



(b): 




Q 2 (GeV 2 ) 



Figure 2: The measured differential cross-sections (a) da/dE ] f? B and (b) da/dQ 2 
for inclusive-jet production (dots) using different jet algorithms. The NLO QCD 
calculations with /xr = E^ B (solid lines) are also shown. The cross sections for 
the anti-k? and kx fS/ algorithms were multiplied by the scale factors indicated in 
brackets to aid visibility. The error bars on the data points are smaller than the 
marker size and are therefore not visible. The shaded bands display the uncertainty 
due to the absolute energy scale of the jets. The lower part of the figures shows 
the hadronisation- correction factor applied to the NLO calculations together with 
its uncertainty (hatched bands) for each jet algorithm; the hadronisation- correction 
factor for the kj- algorithm was shifted by the value indicated in brackets to aid 
visibility. 



23 



o 

-J 

z 

© 



II I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

1-2 :(a ) k T 

« ' * ' 1 f l!!lllilllliniJ;!:llllilli!!llllinillii , 

• ZEUS 82 pb 

I I I I I I I I I I I 



1 [ V/»///?/////M////////J///f/// 

Q_8 — \!i 1 1 1 theoretical uncertainty 

J I I I | I I I I | I I I I | I I I I | I I I I | I I I I i x i i | i i i i 

1.2 7 anti-k T ^^^^^^^^^^ 1 

i = io^¥///^//// / »//////// / i ////w^^///////,v///////,-^ 

0.8 



1.2 \r SIScone 



ZEUS 



miwiiiimiiiiMiiiiiiiw.iiiiiiim 



3 1-2 

? i 

•2 0.8 
" 1.2 
1 

0.8 
1.2 
1 

0.8 




E^ B (GeV) 



(b) k, 

V / T 




1 1 1 1 












mJ 1 — i — iiiii y 

r anti-k T 


1 — Mf 


III -. 


H ; 


E EZSmZZZZZSZzzz^ 

+1 1 1 — Mill h 


1 — Mf 


Ill 1 




r SIScone 






|H ; 


f — ^JSJSMsmMTn 

-.1 i 




,1,1 




^d) anti-k T 




j ; 


-1 1 1 1 Mill 1 


1 

1 1 1 


iii i 


V SIScone "' ' \ 
E mUHJJim^£^-\ -=j - 




if| 1 1 — Mill 1 1 — Mill 1 = 

~- anti-k T - 


=i i i i i = 


2 3 

10 10 




io 4 





io 4 
Q 2 (GeV 2 ) 



Figure 3: The ratios between the measured cross-sections (a) da/dE^ B and (b) 
da/dQ 2 and the NLO QCD calculations (dots). The inner error bars represent 
the statistical uncertainty. The outer error bars show the statistical and systematic 
uncertainties, not associated with the uncertainty of the absolute energy scale of 
the jets, added in quadrature. The hatched bands display the total theoretical un- 
certainty and the shaded bands display the jet-energy scale uncertainty. The ratios 
of the measured cross sections (dots) anti-k^/k^, SIScone/kT and anti-kx /SIS cone 
as functions of (c) E^ B and (d) Q 2 . In these plots, the outer error bars include 
also the uncertainty of the absolute energy scale of the jets. The predicted ratios 
based on calculations which include up to 0{pP s ) terms are also shown (solid lines). 
The hatched bands display the theoretical uncertainty on the ratio. 
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Figure 4: Overview of the theoretical relative uncertainties for the inclusive- 
jet cross sections in the kinematic range of the measurements as functions of (a) 
£y et B and (b) Q 2 for the anti-kx (shaded areas) and SIScone (hatched areas) jet 
algorithms. Shown are the relative uncertainties induced by the terms beyond NLO, 
those on the proton PDFs, that on the value of a s (Mz) and that on the modelling 
of the QCD cascade. 
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Figure 5: The measured differential cross-section da/dE^ for inclusive-jet pro- 
duction in different regions of Q 2 (dots) using the (a) anti-kx and (b) SIScone jet 
algorithms. The measured and predicted cross sections have been multiplied by a 
scale factor as indicated in brackets to aid visibility. Other details as in the captions 
to Figs. H and [3 
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Figure 6: The ratios between the measured differential cross-sections da/dE 1 ^ 
presented in Fig. [5| and the NLO QCD calculations (dots). Other details as in the 
captions to Figs. [H and [3 
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